Signatures of Explosion Models for SN la &: Cosmology 



P. Hoflich 

Dept. of Astronomy, University of Texas, Austin, USA 

Abstract. Based on detailed models for the progenitors, explosions, light 
curves (LCs) and spectra, we discuss signatures of thermonuclear explosions, 
and the implications for cosmology. Consistency is needed to link observables 
and explosion physics. SNe la most probably result from the explosion of a 
degenerate CO- White Dwarf (WD) close to the Chandrasekhar mass. There is 
strong evidence that most of the WD is burned with an extended outer layer of 
explosive C-burning products (0,Ne,Mg) and very little C remaining. Overall, 
the chemical structure is radially stratified. This leads to the currently favored 
delayed detonation model in which a phase of slow nuclear burning as a deflagra- 
tion front is followed by a detonation phase. The importance of pre-conditioning 
became obvious. Within a unified scenario, spherical models allow to understand 
both the homogeneity and basic properties of LCs and spectra, and they allow to 
probe for their diversity which is a key for high precision cosmology by SNe la. 
We emphasize the relation between LC properties and spectra because, for lo- 
cal SNe la, the diversity becomes apparent the combination of spectra and LCs 
whereas, by enlarge, we have to for high-z objects. We show how we can actu- 
ally probe the properties of the progenitor, its environment, and details of the 
explosion physics. We demonstrate the influence of the metallicity Z on the pro- 
genitors, explosion physics and the combined effect on light curves. By enlarge, a 
change of Z causes a shift of along the brightness decline relation because Z shifts 
the balance between 56Ni and non-radioactive isotopes but hardly changes the 
energetics or the ^^Ni distribution. However, the diversity of the progenitors 
produces an intrinsic dispersion in B-V which may pose a problem for reddening 
corrections. We discuss the nature of subluminous SN1999by, and how it can be 
understood in the same framework as 'normal-bright' SNe la. At the example 
of SN200du, we show the influence of the progenitor system and distribution of 
isotopes on light curves. In both objects, we have seen clear evidence for some 
departure from sphericity probably due to circumstellar interaction and stellar 
rotation but the 3D signatures of deflagration fronts remain an elusive feature. 



1. Explosions, Light Curves and Spectra 

The last decade has witnessed an explosive growth of high-quality data for su- 
pernovae. Advances in computational methods provided new insights into the 
physics of the objects, and advances in cosmology. Both trends combined pro- 
vided spectacular results not only for astronomy and the origin of elements but 
also for nuclear, high energy and particle physics, and cosmology. Further im- 
provements and the quest for the nature of the dark energy requires an increased 
accuracy for distance determinations from 10 % to about 2 to 3 % (Weller & 
Albrecht 2001) making evolutionary effects with redshift a main concern, and a 
a better understanding of the physics of SNe la a requirement. There is general 
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Figure 1. Envelopes require radiation hydro dynamical calculations well be- 
yond classical stellar atmospheres or the commonly used 'light bulb' approach 
which assume well defined photospheres with a black body as inner boundary. 
The temperature T, energy deposition due to radioactive decay , Rosseland 
optical depth Tai^ (left scale) and density (right scale) are given as a function 
of distance (in lO^^cm) for a typical SNe la at 15 days after the explosion. For 
comparison, we give the temperature Tgrey for the grey extended atmosphere. 
The two dotted, vertical lines indicate the region of spectra formation. 

agreement that Type la Supernovae (SNe la) are the result of a thermonuclear 
explosion of a degenerate C/0 white dwarf (WD) with a mass close to the 
Chandrasekhar limit. These scenarios allow to reproduce optical/infrared light 
curves (LC) and spectra of SNe la reasonably well. Nowadays, we understand 
the basic, observational features. SNe la appear rather homogeneous because 
nuclear physics determines the structure of the WD, the explosion, light curves 
and spectra: (1) the WD is supported by degenerate electron pressure, (2) the 
total energy production during the explosion is given by the release of thermonu- 
clear energy, and (3) the light curves are powered by the radioactive decay of 
^^Ni produced during the explosion To first order, the outcome hardly depends 
on details of the physics, the scenario, or the progenitor ("stellar amnesia"). 
Homogeneity of SNe la does not (!) imply a unique scenario, and it took the 
revolution in observational methods with respect to time and wavelength cover- 
age to reveal differences and expose the diversity of within SNe la. For recent 
reviews see Branch (1999) and Hoflich et al. 2003). 
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Figure 2. Light curve properties of DD models with various transition 
densities (triangles, Hofiich et al. 2002), and with a WD originating from 
progenitors with main sequence masses Mms between 1.5 to 7 Mq and to 
solar metallicity (open circles, Dominguez et al. 2002). We give the maximum 
brightness My as a function of 56Ni production (upper left) and the tran- 
sition density (lower left), the dependence of the intrinsic color B-V(My) 
on My (upper right), and My as a function of the post-maximum decline 
AMy{15d){V) in V over 15 days (lower left). In addition, we plot observed 
AMy{15d) of 66 well-observed SNe la (Phillips, 2003). For the transformation 
between AMy{15d)(B) to AMy(15d)(V), see Garnavich et al. (2004). 
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Figure 3. Influence of the metallicity Z on the B and V hght curves for a 
progenitor star of 7 Mq on the main sequence based on a delayed detonation 
model. For the progenitor, Z[0/Fe] is taken according to Argast et al. 2000. 




Figure 4. Analysis of the subluminous SN1999by. We show IR-spectra at 
day 11 (upper left), the B and V light curves (right plots), and the chemical 
structure (lower panel). The explosion and evolution of the spectra are calcu- 
lated self-consistently with the only free parameters being the initial structure 
of the exploding White Dwarf, and a parameterized description of the nuclear 
burning front. Without further tuning, the spectra and their evolution with 
time can be reproduced. Up to maximum light, the spectra are formed in 
layers of explosive C-burning, followed by the layers of incomplete Si burn- 
ing. The long duration of these phases provides a lower limit for the initial 
WD mass which is close to the Chandrasekhar mass. Unburned material is 
restricted to the high velocity regime, i.e. the very outer layers. This strongly 
supports DD models for both the normal bright and subluminous SNe la. 

A detailed discussion is well beyond the scope of this contribution. We want 
to illustrate some aspects and, for details, refer to the original publications. 
- Consistency is required to link the observable and the progenitor system, pro- 
genitor and explosion physics. By the physical conditions, consistency is also 
demanded for the treatment of hydrodynamics, rate equations, and radiation 
transport (Fig.l). Density structures require detailed hydrodynamics, low den- 
sities cause strong non-LTE effects throughout the entire envelopes and the radi- 
ation field differs from a black body, chemical profiles are depth dependent, en- 
ergy source and sink terms due to hydrodynamical effects and radioactive decays 
dominate throughout the photon decoupling region, and all physical properties 
are time-dependent because the energy diffusion time scales are comparable to 
the hydrodynamical expansion time scale (Hoflich, 1995). Our approach signif- 
icantly reduces the number of free parameters, namely the initial structure of 
the progenitor, the accretion rate on the WD, and the description of the nuclear 
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Figure 5. High velocity Call feature as tell-tail for interaction within the 
progenitor system. We show the Call IR feature observed in SN 2003du on 
May 6rd in comparison with theoretical models at about 15 days after the 
explosion (upper right), and its evolution with time (left). The calculations 
are based on a delayed detonation model which interacted with a nearby, shells 
of 0.02 Mq and solar composition during the early phase of the explosion. 
The dominant signature of this interaction is the appearance of a persistent, 
secondary, high velocity Ca II feature. Without ongoing interaction, no H or 
He lines are detectable. Note that, even without a shell, a secondary Ca II 
feature can be seen for a period of 2 to 3 days during the phase when Ca HI 
recombines to Ca II emphasizing the importance of a good time coverage for 
the observations. Nearby shells mainly change early time LCs (lower right) 
due to blocking by Thomson optical depth in the shell. In contrast, ongoing 
interaction will change the late time luminosities (from Gerardy et al. 2004). 



burning front. The light curves and spectral evolution follow directly from the 
explosion model without any further tuning. 

- The best current explosion scenario is the delayed detonation (DD) model 
(Khokhlov 1991). These models reproduce reasonable well the color evolution 
of optical and IR-LC, and the spectral evolution of both normal bright (Hoflich 
1995) and subluminous SNe (Fig. 4, Hoflich et al. 2002) within a unified sce- 
nario. By enlarge, DD models produce radially stratified chemical structures for 
elements of explosive carbon, oxygen and Si burning in agreement with the ob- 
servations from the optical and near IR (e.g. Fig. 4, Barbon et al. 1989, Wheeler 
et al. 1998, Marion et al. 2003) The brightness decline relation (Phillips 1993) 
can be understood as an opacity effect (Fig. 3 , Hoflich et al. 1996, Maeda et 
al. 2003) including the color evolution in B-V and the brightness range (Fig. 
2 Hoflich et al. 2002) by varying a single, free parameter for the amount of 
burning when the deflagration turns into a detonation which determines the pre- 
expansion of the WD. In spherical models, this quantity is often parameterized 
by a transition density pt^. Other properties produce a dispersion. E.g. changes 
in the progenitor, i.e. in main sequence mass Mms and Z, alter the explosion 
energy and ^^Ni mass, respectively, but, by enlarge, along the My(^^A^z) and 
AMy(15(i) relations (Hoflich et al. 1998, 2001). However, the diversity in pro- 
genitors cause an intrinsic dispersion in B-V of 0.1^ which, in practice, seriously 
hampers the corrections by interstellar reddening. 
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Figure 6. Late time IR-spectra a probe of the isotopic structure and its 
influence on the optical light curves. On the left, we give the NIR spec- 
trum of the "Branch-normal" SN 2003du on Febr. 27, 2004 with Subaru, in 
comparison with the theoretical line profiles (solid thick line) based on a DD 
model with a central region consisting of non-radioactive iron group elements. 
In addition, the individual components of the forbidden [Fe II] transition at 
1.644 /im are given for the original delayed detonation model (solid) and 
mixed chemistry (light) normalized to the maximum line flux (dotted) and 
the wings (dashed), respectively. Mixing of the inner iron-rich layers of 56Ni 
and stable isotopes (Fig. 4) is to be expected from current 3D models during 
the deflagration phase which is dominated by RT instabilities, and would pro- 
duce round profiles which seem to be at odds with the observations. Possible 
explanations may be that small-scale, pre-existing velocity fields are impor- 
tant for the propagation of nuclear flames. On the right, the visual light curve 
and {B — V) are given for the same delayed detonation model but with and 
without mixing of the inner layers. Differences in V and B-V are ?^0.2^ and 
0.05^, respectively. In effect, mixing redistributes ^^Ni from the outer to 
the inner layers which decreases the photospheric heating at about maximum 
light but increases the 7-trapping later on (from Hoflich et al. 2004). 



- The IR is a primary tool to study intermediate mass elements, and spectropo- 
larimetry and line profiles allow to probe for 3-D signatures (Hoflich 1991, Wang 
et al. 1997, Bowers et al. 1997, Wheeler et al. 1998, Howell et al. 2001, Hoflich 
et al.2002). 

- Our models allow to study particular physics and relations to the properties 
of the progenitors in 'isolation' such as signatures of the progenitor system (Fig. 
5, Gerardy et al. 2004) and the distribution of neutron rich isotopes (Fig. 6, 
Hoflich et al. 2004). In extensive studies, we investigated the influence of the 
central density of the WD, the ignition process, and metallicities (which are 
expected to change with time) (Hoflich et al. 1998, 2001, Dominguez et al. 
2002, see also Timmes et al. 2003). Pre-conditioning of the WD is a key to 
understand the diversity of SNela. Metallicity will effect the statistical sample 
because it changes the life times and radii of stars and, consequently, the binary 
system with mass overflow, and the physics of individual objects. A change of 
metallicity or progenitor MS mass effects both the progenitor evolution, WD 
structure and the observable LC shapes simultaneously, and the combined effect 
must be considered when comparing to observations. As can be seen from Fig. 
3. The change of the initial Fe changes the core He burning during the stellar 
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Figure 7. Analysis of 3D effects in the subluminous SN1999by based on 
flux and polarization data. Upper panel: Comparison of the NIR spectrum 
on May 16 (left) with a spherical, subluminous delayed detonation model 
without and with mixing as expected from detailed 3-D deflagration models 
(Khokhlov, 2001) which does destroy any fit due to excitation of intermediate 
mass elements (S,Si), and the absorption by iron-group elements. Lower, 
left panel: Energy deposition by 7-rays at day 1 (left) and 23 (right) for a 
3D deflagration model based on our full 3-D MC gamma ray transport. At 
about day 23, the energy deposition is not confined to the radioactive ^^Ni 
ruling out dumpiness as a solution to the excessive excitation of S and Si 
lines. Lower, right panel: Optical flux and polarization spectra at day 15 
after the explosion for the subluminous 3-D delayed-detonation model in com- 
parison with the SN1999by at about maximum light. In the observations, the 
polarization angle is constant indicating rotational symmetry of the envelope, 
and an axis ratio A/B of 1.17 (from Hoflich et al. 2002, Hoflich 2002, Howeh 
et al. 2001). 



evolution which determines the C/0 ratio, whereas the '^'^Ne directly effects 
the explosive nuclear burning by shifting the Nuclear Statistical Equilibrium 
slightly away from ^^Ni. The former effect alters the energetics of the explosion 
(explosion energy) whereas the latter changes the energy gain by radioactive 
decay of ^^Ni which powers the LCs. For low metallicities, the influence of the 
energetics dominates the changes in the LC shape. It causes a change Atrise of 
the rise-time for a given decline ratio and, thus, an offset AM in the brightness 
decline relation by AM ^ 0.1 x A{trise)' 

These results do not (!) imply that all processes involved are spherical 
in nature but they are a result of the 'stellar amnesia' mentioned above, i.e. 
the properties depend mostly on integral quantities and basic nuclear physics. 
Despite the successes of spherical delayed detonation models, it becomes in- 
creasingly obvious that multidimensional effects are important towards a better 
understanding of SNe la. 

- Asymmetries in a SN la may be present for a variety of reasons. Within 
the framework of the explosion of a massive WD, reasons and all seems to be 
realized. Reasons include instabilities in the burning front (Khokhlov 1995, 
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2001, Reinecke et al. 2002 Gamezo et al. 2003), namely the deflagration phase, 
off-center detonations or the transitions from deflagration to detonation (Livne 
1999), rapid rotation of the WD, and interaction with the accretion disk and 
the companion star. Whereas we have evidence for the existence of the latter 
effects, direct evidence is still elusive for the RT instabilities or, more precisely, 
such evidence has not been seen although it was expected (Figs. 4 &; 6). 

Both from the remnants and the lack of strong polarization in thermonuclear 
supernovae, limited scattering in brightness, we can conclude that the envelopes 
have rather small global deviations from spherical geometry of the order of 5 
to 20 %. For cosmology, deviations from sphericity will introduce a directional 
dependence of the luminosity of the order of about 4-6 %. 

References 

Argast D., Samland M., Gerhard O. E., Thielemann F.-K. 2000, A&A 356, 873 
Barbon R., Benetti S., Rosiono L., Cappellaro E. Turatto M.1990 A&A 237, 79 
Bowers, E. J. C, Meikle, W. R S., Geballe, T. R., Walton, N. A., Pinto, R A., Dhillon, 

V. S., Howell, , & Harrop-Allin, M. K. 1997, MNRAS, 290, 663 
Branch, D. 1999, ARA&A, 36, 17 

Dominguez, I., Hoflich, R, Straniero, O. 2001, ApJ, 557, 279 
Gamezo, V. N. et al. 2003, Science, 299, 77 
Garnavich et al. 2004, ApJ, in press 

Gerardy, C. L., Hoflich, P., Fesen, R. A., Marion, G. H., Nomoto, K., Quimby, R., Schae- 

fer, B. E., Wang, L., & Wheeler, J. C. 2004, ApJ, in press (astro-ph/0309639) 
Hoflich, P., Gerardy C, Nomoto K., Motohara K., Fesen R.A., Maeda K., Ohkubo T., 

Tominaga N, ApJ, in press (Dec. 20, 2004) & astro-ph/2004 
Hoflich, P., Gerardy, C, Linder, E., & Marion, H. 2003, Lecture Notes in Physics 635, 

Springer Press, p. 203 & astro-ph/0301334 
Hoflich, P., Gerardy, C, Fesen, R., & Sakai, S. 2002, ApJ, 568, 791 
Hoflich, P., Wheeler, J. C, Thielemann, F. K. 1998, ApJ, 495, 617 
Hoflich, P., Khokhlov, A., Wheeler, J. C, Phillips, M. M., Suntzeff, N. B., & Hamuy 

M. 1996, ApJ, 472, L81 
Hoflich, P., Khokhlov, A. 1996, ApJ, 457, 500 
Hoflich, P. 1995, ApJ, 443, 89 

Howell, A., Hoflich, P., Wang, L., Whe eler, J. C. 2001, ApJ, 556, 302 
Kasen D., et al. 2004, ApJ, in press & |astrQ-ph/0301312 
Khokhlov, A. 2001, astro-ph/0008463| 
Khokhlov, A. 1991, ApJ, 245 114 
Livne E. 1999, ApJ 527, L97 

Marion H., Hoflich P., Vacca W.D., Wheeler C.J. 2003, ApJ 591, 316 
Maeda, K., MazzaH, P., Deng, J., et al. 2003, ApJ, 593, 22 
Phillips, M.M., 1993, ApJ, 413, 105 

Phillips, M.M., 2003, Lecture Notes in Physics 635, Springer Press, p. 175 
Reinecke, M., Hillebrandt, W., & Niemeyer, J. C. 2002, AAP, 391, 116 
Timmes F., Brown E., Truran 2003, ApJ, 590L 
Wang, L., Wheeler, J. C, Hoflich, P. 1997, ApJ, 476, 27 

Wang, L., Baade, D., Hoflich, P., Wheeler, J. C, Kawabata, K., Nomoto, K. 2003, ApJ, 
in press 

Wheeler C.J., Hoflich P., Harkness R., Syromillo J. 1998, 496, 908 
Weller J., Albrecht A. 2001, Physical Review Letters, 86, 1939 



